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Abstract 

The 3 P2~ 3 F2 pairing model is generally considered to provide an adequate descrip- 
tion of the superfluid states of neutron matter at densities some 2-3 times that of 
saturated symmetrical nuclear matter. The problem of solving the system of BCS gap 
equations expressing the 3 P2~ 3 F2 model is attacked with the aid of the separation ap- 
proach. This method, developed originally for quantitative study of S-wave pairing in 
the presence of strong short-range repulsions, serves effectively to reduce the coupled, 
singular, nonlinear BCS integral equations to a set of coupled algebraic equations. For 
the first time, sufficient precision becomes accessible to resolve small energy splittings 
between the different pairing states. Adopting a perturbative strategy, we are able 
to identify and characterize the full repertoire of real solutions of the 3 P2~ 3 F2 pair- 
ing model, in the limiting regime of small tensor-coupling strength. The P-F channel 
coupling is seen to lift the striking parametric degeneracies revealed by a earlier sep- 
aration treatment of the pure, uncoupled 3 P2 pairing problem. Remarkably, incisive 
and robust results are obtained solely on the basis of analytic arguments. Unlike the 
traditional Ginzburg-Landau approach, the analysis is not restricted to the immediate 
vicinity of the critical temperature, but is equally reliable at zero temperature. Inter- 
esting connections and contrasts are drawn between triplet pairing in dense neutron 
matter and triplet pairing in liquid 3 He. 



1 Introduction 

It is widely accepted that triplet pairing between constituent spin- 1/2 fermions gives rise to 
superfluid phases in liquid 3 He at millikelvin temperatures and in the neutronic component 
of the quantum fluid interior of a neutron star at temperatures in the hundred keV range 
and below. In both cases, the density is so high that the familiar singlet S-wave gap is 
quenched by the dominant effect of the short-range repulsion in that channel. Instead, 
pairing is favored in a channel with an odd orbital momentum L — 1, and therefore in the 
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triplet spin state. The pairing mechanisms active in these two examples produce interesting 
distinctions between them. In superfluid 3 He, pairing is triggered by spin fluctuations pQ, 
and the B-state (or B-phase) with total angular momentum J = occupies most of the 
superfluid phase diagram. However, the spin-fluctuation mechanism induces relatively tiny 
energy splittings between states with different J values. Consequently, very close to the 
critical temperature T c there exists a phase transition from the B-phase to the A-phase, 
which involves a combination of J = 1 and J = 2 pairing channels. 

The situation in neutron matter is rather different. Referring to the experimental data on 
the energy dependence of the nn scattering phase shifts j2] , one may infer that central forces 
between a neutron pair are quite weak, and that pairing must be due predominantly to the 
spin-orbit force. Based on empirical analyses of spin-orbit splitting in finite nuclei, the latter 
force is expected to be insensitive to polarization and correlation effects. The spin-orbit 
pairing mechanism implies that the spin S and orbital momentum L of the pair cease to be 
conserved separately, and pairing in the J = 2 channel dominates IU 13 EH 13 El E] ■ 

In both Fermi quantum liquids, the effective in-medium interaction contains a component 
that mixes two-body states with orbital angular momenta L ± 2. Specifically, pairing 
channels with L = 1 and L = 3 are coupled by the magnetic dipole force in the case of liquid 
3 He and by the tensor force arising from pion exchange, in neutron matter. However, the 
effect is of minuscule importance in the 3 He problem, since the magnetic dipole component is 
smaller in magnitude than the dominant central part of the interaction by a factor 10~ 7 [TJ. 
By contrast, the magnitude of the dominant spin-orbit force in neutron matter is only a 
few times larger than the tensor component, so the effects of the latter cannot be neglected. 
(More specifically, the parameter measuring the strength of the tensor force relative to the 
spin-orbit force varies around 0.3 in the density interval po < P < 3po [IU E] if in-vacuum 
interaction constants are adopted, p being the saturation density of symmetrical nuclear 
matter.) Moreover, pion exchange is responsible for the most powerful fluctuations in neutron 
matter. Accordingly, it is imperative to give careful attention to 

3 P2^ 3 F2 channel coupling in the quantitative description of pairing in neutron matter. 

Reliable prediction of the phase diagram of triplet pairing over an extensive temperature 
range has proven to be extraordinarily difficult. The traditional tool for elucidation of the 
phase diagram in liquid 3 He and other systems manifesting superfluidity or superconductiv- 
ity has been the Ginzburg-Landau functional approach. But since Ginzburg-Landau theory 
is valid only near the critical temperature, it has little relevance to neutron stars, which cool 
down to temperatures one or two orders of magnitude below T c by a thousand years after 
their birth in a supernova event. There is of course the alternative of brute-force iterative 
solution of the gap equation, which has been widely used in quantitative studies of singlet 
S-wave pairing [TT]. However, this strategy is generally afflicted with slow convergence and 
uncertain accuracy both for S-wave interactions containing strong short-range repulsions and 
for pairing in higher angular momentum states, where one must deal with a multitude of cou- 
pled nonlinear singular integral equations. The limitations of standard iterative approaches 
become particularly serious when one seeks to construct the superfluid phase diagram of the 
system, which is sensitive to tiny energy splittings between the different solutions of the BCS 
pairing problem. 

Explication of the complete superfluid phase diagram of neutron matter has two facets. 
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First, one must identify and characterize the set of admissible solutions of the BCS gap equa- 
tions arising in the 3 P2- 3 F 2 pairing problem, i.e., one must find the "spectrum of solutions." 
Second, one must determine the relative stability of the different solutions under variation 
of density, temperature, and other relevant parameters, so as to uncover the possible phase 
transitions and map out the actual phase diagram. 

The purpose of this paper is to present a detailed account of the substantial progress 
that has been made on the first of these tasks through application of the separation method 
developed in ref. O IE] for robust and accurate solution of BCS gap equations. (The reader 
may consult ref. ^3] for a comprehensive review of the separation approach.) Many of the 
results obtained here have been exposed in condensed form in earlier works (HI E]; it is 
our intent here to offer a more complete justification of these findings, and to provide an 
enhanced understanding of their wider implications. 

In the separation method, the BCS system is recast so as to isolate the major, loga- 
rithmically divergent contributions to the pairing effect and treat them separately from the 
remaining features of the problem, which are very insensitive to the presence of the gap and 
its particular value. This method, which in essence reduces the problem to solution of system 
of algebraic equations, is equally reliable and precise in the limiting regimes T — > T c and 
T — > 0, as well as in between. In concrete calculations, we may accommodate and exploit 
the fact, based on the experimental P-scattering phases [2|, that the central components of 
the in-vacuum nn interaction nearly compensate each other. We assume - quite plausibly 
- that this feature is maintained by the effective nn interaction within neutron matter. We 
furthermore assume - with less confidence - that the smallness of the parameter charac- 
terizing the importance of the tensor force relative to the spin-orbit force is maintained in 
dense neutron matter. It must be acknowledged that medium modification may be more 
significant in the case of the tensor force than it is for the spin-orbit component, especially 
in the vicinity of the phase transition leading to pion condensation ^JEJ]- Investigation of 
this issue calls for a concerted effort and will be pursued elsewhere. 

The current study of the 3 P2- 3 F2 model is exhaustive in its assembly of the set of solutions 
whose structural expression involves only real numbers. The implementation of this program 
requires no numerical computations; remarkably, everything can be done analytically. In 
principle, the same method can be readily applied to the determination of complex solutions. 
However, this more complicated problem inevitably entails some computer work. 

The paper is organized as follows: In Sec. 2, the coupled-channel BCS formalism is stated, 
the separation method is applied to the gap equations, and further analysis establishes the 
explicit equations of an incisive perturbative treatment of the 3 P2- 3 F 2 pairing model. In 
Sec. 3, the parametric degeneracy inherent in the unperturbed 3 P2 problem is lifted, and all 
real solutions of the perturbed model are found. Section 4 furnishes a convenient and succinct 
catalog of these solutions. Finally, Sec. 5 is devoted to a discussion of potential vulnerabilities 
of our approach to the 3 P2- 3 F2 problem in neutron matter, as well as informative connections 
with aspects of triplet pairing in liquid 3 He. 
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2 The Basic Set of Gap Equations 



Fixing the spin and isospin of the pairing state at S = 1 and T = 1 (triplet-triplet), we start 
with the partial-wave decomposition 

A^(p) = E Ai M (p) (1) 

J,L,M 

of the 2x2 gap matrix in terms of the spin-angle matrices 

K») a/3 = E C^f{ti lsM Y LM M) (2) 

Ms Ml 

and the multi-channel BCS gap equations 

UL\J\M\ 

«!«A;*wto (3) 

that embrace the formal problem to be solved [HI El El- The latter equations contain the 
spin trace 

^ JlMl (n) = Tr G J L fi(n)} (4) 

and the interaction matrix elements (p\V^ L ,\pi) appearing in the partial- wave expansion 



V(p, Pl ) = E {- 1 )^WIv\Pi)Glj^) (G#j(ni))* (5) 

LL'JM 



of the block of Feynman diagrams irreducible in the particle-particle channel. The quasipar- 
ticle energy 



E(p) = y/?(p) + D*(p) 

e( P ) + l E (A^(p))*A^(p)^ JlMl (n) (6) 

Z LJMLiJiMi 

is assembled from the single-particle spectrum £ (p) of the normal Fermi liquid, measured rel- 
ative to the chemical potential /x and often parametrized with an effective mass M*, together 
with the gap components /S. J L M (p). We note that E(p) takes on an angular dependence by 
virtue of the spin trace Sj^ lMl (n), which in principle greatly complicates the task of solv- 
ing the system (j3J. However, this angular dependence only comes into play near the Fermi 
surface; hence it can be ignored in those integral contributions to the r.h.s. of Eq. © in 
which the region around the Fermi surface is suppressed [HI El- Consequently, in practice 
the gap equations approximately decouple in the variables L', Li, and Ji, by virtue of the 
orthogonality property 

J S J L ^ (n)dn = 5 LLl 6 JJl 5 MMl . (7) 
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Treatment of the exact, coupled equations © is therefore not as difficult as it appears at 
first sight (though still not at all trivial). 

The angle dependence of the function D(p) makes it awkward to speak of the energy 
gap. Therefore it is conventional to introduce the quantity 

W(k F )} 1/2 (8) 

as a representative measure of the pairing gap in the quasiparticle spectrum, where the 
overbar signifies an angle average. 

As we have argued in the introduction, the spin-orbit component of the neutron-neutron 
interaction exerts a strong influence on pair formation in dense neutron matter, favoring 
the condensation of pairs in the 3 P2 state. Accepting this widely held view, the simplifying 
feature just revealed implies that contributions to triplet pairing from nondiagonal terms 
with L', Li 7^ 1 or J\ ^ 2 on the r.h.s. of Eq. (1) can be evaluated within perturbation 
theory. Defining vp = (pf\Vu\Pf), the relevant coupling parameter is r\ = — (p f \V^\pf)/vf- 
The evaluation is carried out in terms of the set of "principal" gap amplitudes Af M (p), with 
M running from —2 to 2. A further simplification ensues from time-reversal invariance, 
which implies that only three of these five quantities can be independent, say A\ M (p) with 
M = 0,1,2. 

Focusing on the nondiagonal contributions to the r.h.s. of Eq. (1), we observe that 
two of them are of leading significance. The first contains the integral of the product 
V£S^ 2Ml Al M \ and the second contains the integral of the product V&Slff 2Ml A 2 3 Ml . Re- 
stricting attention to these contributions, we arrive at the 3 P2- 3 F 2 pairing problem, which 
has been studied both analytically and numerically in earlier work jHJ El El IS] ■ The list 
of participating states appears to be exhausted upon addition of the 3 Po and 3 Pi pairing 
channels, which are deemed to be of lesser importance; their role will be examined later in 
this paper. 

Adopting the 3 P 2 - 3 F 2 pairing model as circumscribed above, the BCS system (J3J) takes 
the explicit form 

Mi J J 2£ (Pi) 
+ £/ / (p|^bi)^r Ml (n 1 ) tanh ^ 0( r Pl) / 2T) A^Wpldndm , 

Mi J J ^OlPlj 



ai m (p) = e / /(pi^bi)^r M Hn 1 ) tanh ^^ i | /2T) A^( P op 1 ^i^ 1 



(9) 



which will now be subjected to analysis and solution. In writing the r.h.s. of this equation, 
we have replaced the quasiparticle energy E(p; rj) by 

1 /2 

£ (p;r/ = 0) = [e(p)+r%(pT, where D (p) is the gap function of the much-studied 
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3 P2 pairing model in which the tensor coupling between F and P states is ignored. The 
substitution E — > Eq in the integrals on the right is justified by the small relative size of the 
quantities (j>\V? 3 \ Pl ), A 2 3 Ml ( Pl ), and (p\V 3 \\ Pl ). 

A conspicuous feature of the 3 P2 pairing model is the high parametric degeneracy of the 
spectrum of its solutions. This property has been investigated in detail in refs. [HIE] m terms 
of the two ratios Ai = Df^/Df = -D^yfi/D? and A 2 = Df^/Df = D^VQ/D 20 . 
The degeneracy is reflected in the existence of a set of curves Ai(A2) in the plane (Ai,A2), 
upon which all the BCS equations of the 3 P2 problem are satisfied. As we shall see, this 
degeneracy is essentially lifted in the 3 P 2 - 3 F 2 pairing model where rj ^ 0. A finite set of 
points (Ai, A 2 ), depending somehow on the 77 value, replaces the set of solution curves Ai(A 2 ) 
of the 3 P 2 model. 

It is the objective of this article to identify the different solutions of the system @ and to 
establish their structure in the realistic case of small 77. The analysis is aided by the fact that 
the parameters Ai = fi{f]) and A 2 = / 2 (^) are continuous functions of the coupling constant 
77. This property implies that the number of solutions of the 3 P 2 - 3 F 2 pairing problem as well 
as their structure remains the same no matter how small 77 is. Consequently, implementation 
of our program reduces to determination of the functions /1 and / 2 at 77 = 0, i.e., Ai(t7 = 0) 
and A 2 (?7 = 0). It then becomes apparent that the quasiparticle energy E{rj) may be replaced 
by E on the left-hand-sides of Eqs. as well as on the right, since taking into account 
the difference E{rj) and Eq within (JUJ) cannot, by itself, lift the parametric degeneracy. This 
conclusion is confirmed in the numerical calculations we have performed. 

The nondiagonal integrals on the r.h.s. of Eqs. are rapidly convergent, with the over- 
whelming contributions coming from momenta adjacent to the Fermi surface. This feature 
greatly expedites application of the perturbation strategy. For E(p) significantly in excess of 
the energy gap A^ of Eq. (JBJ), the energies E(p) and |£(p)| are coincident to high precision, 
such that the angular integration in Eq. (jHJ) yields a null result. Thus, when treating the 
nondiagonal contributions it is sufficient to know the minor gap components A 2M (p) at the 
point p = pf, which may be efficiently evaluated in terms of the coefficients D\ M = A\ M (p F ) 
(with M = , 1 , 2). In this process, we retain, on the r.h.s. of the last of Eqs. ©, only the 
dominant contribution containing a large logarithmic factor L = ln(e^/Aj?), where e F is the 
Fermi energy. This factor is angle-independent; therefore the respective angular integral is 
freely evaluated, giving rise to the simple connection 

Af \p = p F ) = -L(p f \VUpf)D\ m = V v F LDl M ~ 7]D\ M . (10) 

In obtaining this relation we have employed the equality 1 = v F L, which holds when one 
keeps only logarithmic contributions. Analogous linear relations are obtained for the other 
minor components A J L M of the gap function (notably A® and A} M ). 

Insertion of the result (fTU|) into the first of Eqs. © leads to the closed system of equations 

a?» + E / / (p!^ 2 i bi)^T 2Ml W tanh ffiM /20 A 2M, {pi)pldpidni 



Mi 



Mi J J 
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X 



for finding the set of three gap functions A\ M (p) with M — 0, 1, 2. 

At this point we invoke the separation method and assert the decomposition H] 



2M„ 



(12) 



of the gap component into a "universal" shape factor x{p) that is independent of the magnetic 
quantum number M and a numerical coefficient that embodies the dependence on M. 
The function x{p)i normalized by x{Pf) = 1, is the solution of a linear integral equation. 
As argued in refs. |Hl E] , this decomposition holds to high accuracy in the problem domain 
under consideration. Accordingly, our problem reduces to the determination of the three 
key coefficients Df M , which obey a set of coupled algebraic equations obtained by setting 
p = pf in Eqs. flTTj) . With M = 0, 1, and 2, these equations read 



Dl M + v F Y,Df Ml 

Mi 
M, 



Vi2M2Mi/ \ | q2M2Mi, 
D 31 \ n ) + D 13 1 



n 



-K"o( n )^ n 



(13) 



with (p{p) = (pIV^IpfV^f, and 



^o(n) 



(Pl^ilbi) 



tanh( J B (p)/2T) 
2^o (p) 



X(p)p 2 dp. 



(14) 



It should be remarked that the integral (J 14)) contains a constant coming from regions lying 
far from the Fermi surface, but this constant does not contribute to the angular integration 
on the r.h.s. of Eq. (fT3*|) . 

In the present work, the search for solutions will be confined to those with real coefficients 
Df . The structure of the phases having complex coefficients Df M can be explored and 
established along the same lines, although the calculations become much more cumbersome. 
To streamline the task of finding solutions, it is helpful to rewrite Eqs. ()13|) in terms of the 
ratios Ai and A2. This step ensures coincidence between the left-hand sides of these equations 
and their counterparts in the model of pure 3 P2 pairing solved in ref. [H]. Substitution of the 
explicit form of Sf;f /2Ml (n) into Eqs. (JT3J), followed by straightforward algebra, gives a system 
of three equations 

A 2 + vf [A 2 (^o + ^5) - A1J1 - J 3 ] = r)v F r 2 , 
Ai + vjr[-(A 2 + l)Ji + Ai(Jo + 4J 5 + 2J 3 )/4] = r]v F r x , 
1 + v F [-(A2J3 + AiJi)/3 + J 5 ] = r]v F r , 



for these ratios and the gap value A F . Here, as before jSj, 

Ji = 



tanhE (p)/2T p 2 dpdn 

fi{#, <p)<p(p) — / x — xipy 



2E (p) 



4tt 



(15) 



(16) 
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with f = 1 - 3z 2 , ft = 3xz/2, f 3 = 3(2x 2 + z 2 - l)/2, and 

f 5 = (1 + 3z 2 )/2 and z = cost?, x = sin?? cos and y = sin^sinyj. Among the integrals Jj 
(i = 1, • ■ ■ 5), only J 5 contains a singular principal term going like \n(ej?/ Ap). All the other 

1 /2 

Ji converge close to the Fermi surface, where E (p) = [£ 2 (p) + -D 2 (n)] with 

k2 



£>o>) 



A 2 

F 



2[l + (A? + A|)/3] 



A 

1 + 3 cos 2 + A 2 , sin 2 + -f (1 + cos 2 . 



1 



2Ai(l + A 2 ) cos 9 sin cosy? + - (A 2 — 4A 2 ) sin 2 9 cos! 



In terms of the coefficients Ai, A 2 , the right-hand sides of Eqs. (fTKJ) read 

r 2 = A 2 s 22 + Ais 2 i + v6s 2 o , 
n = A 2 s i2 + AiSn + V6s 10 



where 



SM2 
SMI 
%0 



2M2,-2, 



n 



^T 21 (n) + 5 1 2 f 21 (n)-5 31 



2M2,-1 



(n) - 5; 



2M2,-1, 



13 



n; 



K (n)dn 
K (n)dn 



5 3 2 f 20 (n) + S 2M20 < 



'13 



n 



-K"o( n )^ n • 



(17) 



(19) 



3 Real Solutions of the 3 P 2 - 3 F 2 Problem 

Equations (fT5j) have three familiar one-component solutions with definite magnetic quantum 
numbers M — 0,1, and 2. To uncover the structure and the spectrum of the multicomponent, 
mixed-M solutions of the perturbed problem, a two-step transformation is applied to the 
system (JT3J). The integral J 5 introduces the gap value into the description, but it is 
irrelevant to the phase structure. As a first step, we combine Eqs. (J15|) so as to eliminate 
terms involving J 5 from the first pair and, at the same time, reduce the number of the J« 
integrals in each equation to two. The resulting equations are 

(A 2 + 1)[3A 1 (A 2 + 1)J -2(A 2 -2A 2 + 6)J 1 ] = rjBx , 

(A 2 + 1)[(A 2 -4A 2 )J 1 + A 1 (A 2 + 1)J 3 ] = VB 2 , (20) 

with 

B 1 = 2A 1 (2A 2 + 3)r 2 -4(A 2 -3)r 1 -6A 1 (A 2 + 2)r , 

B 2 = -Air 2 + 4A 2 r x - 3AiA 2 r . (21) 

Two cases are to be distinguished. The first corresponds to A 2 = —1, a particular 
solution of the pure 3 P 2 pairing problem [S]. In this case, the left-hand sides of Eqs. (J2*U|) 
vanish identically, and hence so must their right-hand sides, leading to the single restriction 

Air 2 (Ai; -1) + 4n(A i; -1) - 3Air (Ai; -1) = . (22) 
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As will be seen, this condition is satisfied at any Ai. Therefore the particular solution 
A2 = — 1 found in the pure 3 P2 problem survives intact when the 3 F 2 coupling is switched 
on. 

Let us now assume A 2 7^ — 1 and proceed to the second step. Following ref. [SJ, we 
perform a rotation 

J z = tcos(3 + usm/3 , . , 

I x = — t sin (3 + u cos (3 , 

with the objective of removing the integral Ji from Eqs. (J2Uj) . To this end, the parameter 
£ = tan (3 is taken as a root of the quadratic equation 

C 2 + ^k-l = 0. (24) 

The integrals Jj transform according to 

Jo — * J (cos 2 /3 — - sin 2 /3) — J 3 sin 2 /3 , 



3 1 
J 1 — > — ^ Jo sin 2 + J3 (cos 2 /3 + - sin^ 

3 1 

Ji -> (-J + -J 3 ) sin/3 cos/?. (25) 



^4 u 2 

Substitution of the transformed integrals into Eqs. (j2U|) yields 



(A 2 + l)[AiJ + A 2 J 3 ] = rjBi, 

(X 2 + 1)[A 1 J + A 2 J 3 ] = -2 V B 2 , (26) 



with 



Ax = -A 1 (1 + A 2 )(2-C 2 )--(A 2 -2A 2 + 6)C, 

A 2 = -3A 1 (1 + A 2 )C 2 -(A 2 -2A 2 + 6)C. (27) 

(Some details of this step are provided in refs. [El EH-) 

The left-hand members of the two equations (j2*o^) are seen to be identical. In fact, the 
universalities of pure 3 P 2 pairing revealed in ref. [S] stem from this key property. The 
solutions of the restricted problem derived from Eqs. ()26|) at 77 = fall into two groups 
composed of states that are degenerate in energy. This remarkable feature is independent 
of temperature, density, and details of the in-medium interaction. There is an upper (i.e., 
higher-energy) group consisting of states whose angle-dependent order parameters have nodes 
and a lower group without nodes (cf. ref. [TS]). In addition to the energy degeneracies, the 
multicomponent pairing solutions, which obey the relation 

(A 2 - 2A 2 - 6A 2 ) (A 2 + 2 - 2A 2 ) = , (28) 

manifest a parametric degeneracy with respect to the coefficient ratios Ai and A 2 , as they in 
general define curves rather than points in the (Ai, A 2 ) plane. 
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Among solutions of the 3 P2- 3 F2 pairing model there exist some for which B\ and B 2 in 
Eq. (J2fijl vanish simultaneously at a certain set of parameters Ai and A 2 satisfying Eq. (|25jl. 
Indeed, consider the set of values 

Ai = 0; A 2 = ±l,±3, (29) 

which determine two-component solutions of the pure 3 P2 pairing problem t 8j. In these cases, 
it may be observed from Eqs. (J21j) and (|18|) that 

#i(Ai = 0, A 2 ) ~ B 2 (X 1 = 0, A 2 ) ~ r x (Ai = 0, A 2 ) 
= A 2 s 12 (A! = 0,A 2 ) + v / 6s 10 (A 1 = 0,A 2 ), (30) 

where Si 2 (Ai = 0, A 2 ) and Si (Ai = 0, A 2 ) are defined by Eq. (|T9"j) . As seen from Eq. (fTTj) . 
the quantities -D (Ai = 0, A 2 ; n) and -Ko(Ai = 0, A 2 ; n) are even functions of cosy? El, 
while both S 3 T 2 (n) + S 2 3 122 (n) + S 31 2 '- 2 (n) + S 2 3 12 '- 2 (n) and S 3 T°(n) + S 2 3 120 (n) appear 
to be linear in cos<^ (for details, refer to the appendix). As a result, both the matrix 
elements si 2 (Ai = 0, A 2 ) and sio(Ai = 0, A 2 ) vanish identically when the integration over ip 
is performed. Hence both of the quantities B\ and B 2 turn out to be zero, and we find that 
Eq. (31) defines a discrete set of valid two-component solutions for the 
3 P 2 - 3 F 2 pairing problem as well as for the pure 3 P 2 case. No other two-component solutions 
of the former problem have been found in numerical calculations. 

Continuing our exploration, suppose that B\ and B 2 do not vanish. Equating the right- 
hand members of Eqs. ()26|) in the 77 — limit, we then obtain an additional relation between 
the parameters Ai (77 = 0) and A 2 (?7 = 0) similar to (|2"2"j). namely 

A 1 r 2 (A 1 , A 2 ) - (A 2 - 3)r 1 (A 1 , A 2 ) - 3A 1 r (A 1 , A 2 ) = , (31) 

the quantities r M being defined by Eq. (fTSj) . Inserting the explicit expressions for the r M , 
this auxiliary condition can be recast as 

G(A 1; A 2 ) 

K (x, y, z; A 1; A 2 )^(x, y, z; X 2 )5(l - x 2 - y 2 - z 2 ) ^ * = , (32) 

271 

where ^(x, z; Ai, A 2 ) = XiR 2 (x, z) — (A 2 — 3)Ri(x, z) — 3Ai-R ( a; 5 z ) an d the quantities Rk are 
given in the appendix. Substituting for the Rk, we obtain 

#(x,z;Ai,A 2 ) 

= 14Ai [5(A 2 + 3)z 4 - 10\ 2 x 4 + 15(A 2 - 3)x 2 z 2 + 6(A 2 + l)(a; 2 - z 2 ) 
+ 70 [(A 2 -3)(A 2 + 3) + 2A 2 " 



xz 3 



+ 70 



2A 2 (A 2 - 3) - 2\j\ x A z - 96(A 2 - 3)(A 2 + l)xz . (33) 

The relation (J3~2"j) supplements the spectral condition (|2~Hj) . As a direct consequence, the 
strong parametric degeneracy inherent in pure 3 P 2 pairing is lifted in the case of 3 P 2 - 3 F 2 
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pairing. With the exception of the straight-line solution A2 = 1 noted above, the solutions 
of the problem are now represented by a set of isolated points in the (Aj, A2) plane. 

The system formed by Eqs. (|2"5j) and (j31|) is amenable to analytic solution. We begin the 
search for solutions of Eq. (J32)) with the particular solution (Ai, A2 = —1) of the pure 3 P2 
pairing problem. In this case, 



4 1 4 r 2 2 
~ Z + X — OX Z 



3(1 



2\2 



while, as seen from (fTTj). the gap function -Dq(A2 
Dl(x,y,z;X 1 ,X 2 = -1) 

Since the integrals 



y-)--A(l-y 2 )(2z 2 -l + y 2 ), (34) 
•1; n) depends only on the single variable 



5(i 

2 V 



y 



jz 4 -3(l -y 
y/1 — y 2 — z' 2 



2\2 



dz 



and 



2z 2 - 1 + y 2 
y/1 — y 2 — z- 



dz 



(35) 



(36) 



both vanish, the integral G(Ai, A 2 = —1) also vanishes at any Ai. Accordingly, the degenerate 
solution A2 = — 1 of the uncoupled 3 P2 pairing problem survives when the 3 F 2 channel is 
involved. 

In search of other solutions, we again make use of the transformation ([23)1 . applying it 
now to the whole integrand of Eq. (J32|) . As found in ref . [Hj , the gap function D reduces to a 
function of the single variable t under such a transformation; hence the same property holds 
for the factor Kq in Eq. fl32(l . which is defined as a functional of Dq by Eq. ()14j) . To ascertain 
how the factor \l/ is transformed, let us write down the results of the transformation for 
simple terms entering this function. In implementing the transformation we omit odd-power 
terms u and w 3 , which do not contribute to the integral ()32|). One obtains 



x 4 - 


-> t 4 sin 4 p + u 4 cos 4 (3 + 6t V sin 2 cos 2 (3 , 






z 4 - 


-> t 4 cos 4 (3 + u 4 sin 4 f3 + 6tV sin 2 (3 cos 2 /? , 






x 2 z 2 - 


-> t 4 sin 2 /3 cos 2 /3 + u 4 sin 2 /5 cos 2 (3 + t 2 u 2 (sin 4 /3 + cos 4 


/?- 


4 sin 2 /3 cos 2 (3) , 


xz 3 - 


-> m 4 sin 3 /3 cos /3 - t 4 cos 3 /? sin /3 + 3t 2 w 2 (cos 3 /3 sin (3 - 


sin 3 


/3 cos /3) , 


x 3 z - 


-> u 4 cos 3 f3 sin /9 - t 4 sin 3 /3 cos /3 - 3t 2 u 2 (cos 3 /3 sin /3 - 


sin 3 


(3 cos /?) , 


2 2 

X — z - 


-> (t 2 -u 2 )(cos 2 /?-sin 2 /3). 







After inserting these relations into the formula (J33|) for \]/, simple algebra yields 



(t, u; Ai, A 2 ) = (1 + CT 2 t^ 4 + Tt 4 + V> 2 - i 2 ) + W^ 2 t 2 



where 



C/(Ai,A 2 ) 



(37) 



(3* 



70|A!(A 2 + 3)C 4 + [(A 2 - 3)(A 2 + 3) + 2A 2 J C 3 
+ 3Ai(A 2 - 3)C 2 + f2A 2 (A 2 - 3) - 2\\\ ( - 2AiA 2 



(39) 
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T(Ai,A 2 ) = 70|-2A 1 A 2 C 4 -[2A 2 (A 2 -3)-2A?]C 3 + 3A 1 (A 2 
- [(A 2 - 3)(A 2 + 3) + 2A 2 ]C + Ai(A 2 + 3)} , 



V(X 1 ,X 2 ) = 24 



7A 1 (A 2 + 1)(1 + C 2 )(1-C 2 ) 
8(A 2 -3)(A 2 + 1)(1 + C 2 )C 



W(X 1 ,X 2 ) = 210|2A 1 (A 2 + 3)C 2 -4A 1 A 2 C 2 + A 1 (A 2 -3)(C 4 + 
jA 2 -3)(A 2 + 3) + 2A 2 ] (C-C 3 ) 
2A 2 (A 2 -3)-2A 2 l (C 3 -C)|. 



+ 
+ 



These results can be simplified slightly by employing the connection 
Ai(C 2 — 1) = (A 2 — 3)C, and one finally arrives at 



U(X 1 ,X 2 ) = 70 



Ai(A 2 + 3)C 4 + 5Ai(A 2 - 3)C 2 - 2AxA 2 



+ (A 2 - 3)(A 2 + 3)C 3 + 2A 2 (A 2 - 3)C 

T(Ai, A 2 ) = 70 Ai(A 2 + 3) + 5A!(A 2 - 3)C 2 - 2A X A 2 C 4 
- (A 2 - 3)(A 2 + 3)C - 2A 2 (A 2 - 3)^ 

V(X 1 ,X 2 ) = 180A!(A 2 + 1)(1-C 4 ), 
W(X 1: X 2 ) = 420A!(A 2 -3)(C 4 -6C 2 + 1). 



The auxiliary integrals 



/ 


-t-' 




du 




-t 2 


Vi 


-t 2 -u 2 


/ 


-t-' 




u 2 du 


VI- 


-t 2 




-t 2 -u 2 



n(l-t 2 ) 



and 



- I 



u 4 du 
VI - t 2 - u 2 



37r(l-t 2 ) 2 
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are helpful in completing the evaluation of G(X±, A 2 ) by integration over the new variables. 
The integrals I , J 2 , and J 4 are related by 

h = |(1 ~ t 2 )h = |(1 - t 2 ) 2 /o , h = \{l ~ t 2 )I ■ (48) 

Using these connections, one can easily verify that integration of the combinations 8m 4 — 
3(1 — t 2 ) 2 and 2m 2 — 1 + t 2 over u gives zero. It follows then that if we make the following 
replacements 

u 2 ^l(l~t% „*->|(l-t»)' (49) 
in Eq. the new function 

*'(Ui,A 2 ) = ± c2)2 { [3U(X 1 , A 2 ) + 8T(A 1; A 2 ) - 4W(X 1 , A 2 )] t 4 

- [6U{X 1 , A 2 ) + 12y(A l5 A 2 ) - AW{Xt, A 2 )] t 2 

+ 3U(X ll X 2 )+W(X l ,X 2 )} (50) 

will guarantee the same result as given by ^ upon integration of Eq. (J32J) . For this integral 
to vanish identically, with the function Ko(t) regarded as arbitrary, the coefficients of all 
powers of t in the function ty' must be zero: 

3U(Xx, A a ) + 8T(A 1; A 2 ) - AW(X X , A 2 ) = , 

3U{X U A 2 ) + 6V r (A 1 , A 2 ) - 2W r (A 1) A 2 ) = , (51) 
3f/(A 1 ,A 2 )+4V r (A 1 ,A 2 ) = 0. 

This system reduces to the chain of equations 

T(A 1; A 2 ) = V(Xx, A 2 ) = W(X U A 2 ) = -^(A 1; A 2 ) , (52) 



where Ai and A 2 are to satisfy the relation ()28)1 . It can be proved that to determine all 
solutions of this set it is sufficient to solve the equation 

V(Ai, A 2 ) = W(Xi, A 2 ) , (53) 

which has the explicit form 

3(A 2 + 1)(1 - C 4 ) = 7(A 2 - 3)(C 4 - 6C 2 + 1) , (54) 



and then to verify that the other equalities in (j52j) hold for these solutions. 
To solve Eq. (|54j) , consider the first branch of Eq. (|28|) , 

A 2 = 2A 2 + 6A 2 , (55) 

for which 

( = tan/5 =^1. (56) 
Ai 
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Then Eq. is recast to 

(A 2 -3)(5A 2 ; + 24A 2 -9) = 0. (57) 

One root of this equation is obviously A 2 = 3, which corresponds to Ai = 6. Another pair of 
roots is given by 

A 2 = - (±V21-4) « 0.350,-5.150, (58) 

yielding respectively 

Ai = ^ V / 2( 17 T3V21) « 1.530, 4.705 . (59) 
Now consider the second branch of Eq. ([2*51 . 

\2 



with 



In this case, Eq. (j^U) becomes 



M = i + 1, 




(A 2 -3)(A^-26A 2 + 29) = 0. 



(60) 
(61) 

(62) 



Its roots are (A 2 , A 2 ) = (3, 2) and 

A 2 = 13 =F 2^ « 1.168, 24.83, 

\ 1 = 2\J 6 =F "s/35 pa 0.579, 6.904 . (63) 
Insertion of these solutions into other ()52|) results in a chain of identities, as required. 

4 Phases of 3 P 2 - 3 F 2 Pairing 

The set of solutions of the 3 P 2 - 3 F 2 pairing problem revealed by the above analysis is depicted 
in Fig. 1 and cataloged in Table 1. 

Figure 1 refers to the multicomponent solutions and represents them by their coordinates 
(Ai, A 2 ) in the two-dimensional parameter space. Only the right half of the Ai — A 2 plane is 
plotted, because the pairing energies are independent of the sign of Ai. The most important 
message of this figure is that (with the exception previously noted) the solution curves that 
represent parametrically degenerate solutions of the 3 P 2 problem shrink to discrete points 
as the degeneracy is lifted by the perturbation that admixes the 3 F 2 channel. 
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1 1 
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Ai 

Figure 1: Parameter sets (Ai,A2) defining the multicomponent solutions of the 3 P2- 3 F2 
pairing problem. Solution curves or points of the pure, uncoupled 3 P2 problem are identified 
with an open circle [or a filled circle] according as their order parameters are nodeless [or 
display nodes]. 



The solutions (which correspond to pairing states and ultimately to superfluid phases) 
divide into two categories: those whose order parameters contain nodes and those whose 
order parameters are nodeless. It is convenient to identify the nodal states with the symbol 
X and the nodeless states with the symbol O. 

First there are the three well-known "one-component" solutions, belonging to magnetic 
quantum numbers M = ±2, M — ±1, and M = 0, respectively. In addition, we have 
established the existence of ten multicomponent solutions that mix states with different 
values of |M|. These are comprised of five nodeless solutions and five exhibiting nodes. 

The five nodeless multicomponent solutions include: 

(i) Two two-component solutions denoted 0±3, which are identical to the "particular" 
solutions found in the pure 3 P 2 case having Ai = and A2 = ±3. 

(ii) Three three-component solutions associated with the branch (|53|l of Eq. (|28j) . Two of 
these, denoted Oi and O4, derive respectively from the upper root of the pair (JH%|) and 
the root A2 = 3. The third, named O2, derives from the lower root of the pair (|58|). 

The five nodal solutions consist of: 

(i) Two two-component solutions X±i that are identical to the "particular" solutions found 
in the pure 3 P2 problem having Ai = and A2 = ±1. 

(ii) Three three-component solutions associated with the branch (|60|) of Eq. (|28|). The 
solutions X 2 and X 4 derive respectively from the lower and upper roots of the pair 
(jHSJ), while X3 derives from the root A2 = 3. 
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Table 1: Identification of the thirteen solutions (or superfluid phases) of the 3 P2 pairing 
model, in terms of the parameters Ai and A2 defining their magnetic content, and in terms 
of their nodeless (states labeled O) or nodal (states labeled X) character. 



Phase 


Ai 


A 2 




M = 


-\r 


M = 1 




M = 2 







3 







-3 


x x 





1 


X_! 





-1 


Ox 


§^2(17-3^21) 


f(V2l-4) 


o 2 


1^/2(17+3^21) 


-§(^21+4) 


4 


6 


3 


x 2 




13-2 v / 35 


x 3 


2 


3 


x 4 


2^/6+^ 


13+2^ 



Appealing to continuity in the parameter 7], it may be expected that the general features 
of the spectrum of solutions delineated by the perturbative analysis will persist even when 
i] is not especially small. 

5 Discussion 

In the present article, we have employed the separation method of ref. jH] to derive, from the 
system of nonlinear BCS integral equations for the components A^ M (p) of the gap function, 
a set of nonlinear algebraic equations for the coefficients specifying the angular structure of 
the superfluid states of interest. Based on this stratagem, we have been able to find all the 
real solutions of the 3 P 2 - 3 .F 2 pairing problem in the regime of vanishingly small rj. Their 
salient feature is a remarkable independence of the temperature T and of any details of the 
interaction V. In principle, the full collection of complex solutions can also be obtained along 
the same lines, but the calculations become much more cumbersome. 

It must be noted that the status of the 3 P 2 - 3 F 2 pairing model, in which contributions 
from 3 P 2 ^ 3 Po or 3 P 2 ^ 3 P! transitions are assumed to be unimportant, is somewhat vul- 
nerable. This assumption should be a safe one in the low-density regions of a neutron star. 
However, it may not hold in the denser core region, where (i) the cancellation between dif- 
ferent contributions to the 

P-wave central component of the scattering amplitude becomes questionable and (ii) the 
amplification of the tensor force due to the pion-exchange renormalization becomes pro- 
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nounced. These effects, which increase with density, may overwhelm the spin-orbit compo- 
nent in the effective pairing interaction, whose strength depends only mildly on density. The 
3 Po pairing channel would then take command and give rise to a superffuid state analogous 
to the B-state of liquid 3 He. In the less challenging situation where 3 P2 — > 3 Po transitions 
can be described within perturbation theory, their effects reduce to a renormalization of the 
matrix elements smm 1; which acquire an additional term of the form 

Ss MMl ~ J S 1 2 f 00 (n)ir (n)rfn J S^(n)K (n)dn . (64) 

An analogous modification occurs when contributions from 3 P2 — > 3 Pi transitions are taken 
into account perturbatively. 

In short, the introduction of new transitions among two-body states results in the mixing 
of components having different values of total angular momentum J. This effect depends 
very strongly on the proximity to the critical temperature T c . It is well known that in the 
immediate vicinity of T c , all the BCS equations decouple; hence the various phases that 
appear may be classified according to the quantum number J. These considerations apply 
as well to superffuid 3 He. The magnetic dipole interaction between the 3 He atoms - which 
is responsible for such a phase separation close to T c - triggers the planar state with J = 1 
[H EOj • This result is unaffected by strong-coupling corrections. 

Let us now turn to a more detailed analysis of the impact of transitions between pairing 
states. We start with the two-component solutions (ffiJj) . namely Ai = 0, A2 = ±1,±3, and 
recall that in the 3 P 2 - 3 P 2 model, both of the quantities Pi(Ai = 0) and B 2 (Xi = 0), which 
behave as ri(Ai = 0), vanish identically. As will now be argued, this property holds in an 
extended version that allows for transitions between states with different J values. Upon 
setting Ai = in the equations (J21)) and (|18|) determining B\ and B 2 , there apparently 
remains a single term containing r^Ai = 0, A 2 ) = A 2 <5si 2 (Ai = 0, A 2 ) + v / 65si (Ai = 0, A 2 ). 
But upon closer inspection of the integrands in this term is also found to vanish. 

According to Eq. (JHJ), -Ko(Ai = 0, n) is an even function of cosy? independently of A 2 , 
while S^ 00 ^), the explicit form of which is given in the appendix, is an odd function of 
cos (p. Nullification of the respective integrals therefore occurs upon integration over (p. 
Consequently, the equality Bi = B 2 = applies even when the transitions 3 P 2 — - > 3 Pi are 
incorporated into the description. It follows that the spectrum (J29|) of the two-component 
solutions of the 3 P 2 - 3 F 2 pairing model endures upon inclusion of the 3 P 2 ^ 3 Po and 3 P 2 — > 3 Pi 
transitions. 

The opposite situation is encountered when one considers the particular solution of the 
3 P 2 - 3 F 2 pairing model defined by the straight line A 2 = —1. The single point Ai = 0, A 2 = — 1 
of this line does survive as a solution of the pairing problem, but the other points do not. 
To confirm this assertion, we first observe that in Eq. (jUJ), contributions from the corrected 
values of r 2 and tq are proportional to Ai and hence vanish provided Ai = 0; further, as 
seen from Eq. (j!8|) . the corrected value of r^Ai = 0) is zero, since contributions from the 
terms of Eq. (|64j) vanish at Ai = 0. To reveal other viable points on the straight line, 
one looks for the zeroes of the terms in the l.h.s. of Eq. (J22|) introduced by 3 P 2 — > 3 Pq 
and 3 P 2 — > 3 Pi transitions. Direct calculation shows that the term corresponding to the 
3 P 2 — > 3 Pi transition vanishes identically on the line A 2 = —1, while that for the 3 P 2 — > 3 Pq 
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transition is proportional to Ai Jg, where the integral Jq is given by Eq. f)16|) . It can be 
checked that this expression has no new zeroes on the straight line A2 = — 1 in addition to 
the point Ai = 0. 

Close inspection reveals that the revised locations of the other solutions found above 
depend on the matrix elements of the interaction V; if the ratio rjp/r) is small then the shifts 
from the old positions are also small. Naturally it is of interest to trace the trajectories 
of solutions as this ratio is varied, thereby exploring the phenomenon of triplet pairing in 
the continuum of cases from that of dense neutron matter to that of superfluid 3 He. This 
investigation will be a subject of future work. 

Our method is to be compared with Ginzburg-Landau (GL) theory, which is generally 
regarded as the standard technique for mapping the spectrum of the phases of systems with 
triplet pairing. In the GL method, the search for diverse phases is based on the construction 
of a suitable free energy functional up to terms of fourth (or even sixth) power in the gap 
value A. This approach allows one to simultaneously evaluate the splitting between the 
different phases and efficiently determine the phase diagram. Another advantage of the 
GL procedure resides in the facility of including strong-coupling corrections ^H! arising 
from the dependence of the effective interaction V on the gap value, an effect that becomes 
important close to the critical temperature T c . Unfortunately, the GL method fails when 
the temperature T is significantly different from T c . Its application to the phase-spectrum 
problem makes a sense only if the phase structure of superfluid neutron matter is independent 
of T. However, special conditions must be met for this to be true in the face of the explicit 
appearance of the factor tanh(.E/2T) in the set of BCS equations, and the explication of 
these conditions is impossible within the GL approach itself. Our method is free of these 
shortcomings. It is equally reliable close to T c and at T = 0. Moreover, the incorporation 
of strong-coupling corrections reduces to the insertion of new terms in V that depend on A 
itself; if this dependence is appropriately specified, no further hurdles must be overcome to 
fully elucidate the triplet superfluid phase diagram. 

In a sequel to this paper, we shall report the findings of a quantitative treatment of the 
transitions between the phases arising in the 3 P2- 3 F2 model and construct the corresponding 
superfluid phase diagram of dense neutron matter. This task requires significant numerical 
effort beyond the analytic developments of the present work. 
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A Explicit Formulas for Matrix Elements 

To find analytic solutions of the system (|5|L explicit formulas for the set of matrix elements 



S 



JMJiMx 
LLi 



n 



1 



x 



;-l) M+1+(L - il+ ^ )/2 — J2 V( 2K + !)( 2J + !)( 2J i + l )( 2L + !)( 2L i + !) 



(-l) (K+|K|)/2 P^ K (^,^)e-^. (65) 





L 




o 





o0( 





J 1 




Ji 




K Mi 


—M J 




L 





are needed, where the Vkk^, 1 ^) are associated Legendre polynomials. Lengthy algebra 
involving 3j- and 6j-symbols yields 



(sT + S 

(sT + s 
(s^- 1 + s. 
(sT + S: 

( o202,-2 , q202, 
I D 31 

(sT + Si 

+ S 
[ST' 1 + S; 

(sT + s 



2020 
31 



2021 
31 



202,-1 
31 



2022 
31 



2120 
31 



2121 | q2121 
13 "T °31 



212,-1 
31 



2122 
31 



/q212,-2 q212,-2 
\&13 + ^31 



2220 
31 



2221 
31 



(sT + s 
(sT + s 

/r<222,-l Q 222,-l 
^D 13 -+- D 31 

(sT + Si 

(^222, — 2 | r-*222, 
^13 + *->31 



2222 
31 



— f 12 sin $ cos sin d cos 3 flA e Jip , 

7vr V 4 / 

— — (l2 sin $ cos $ — sin $ cos 3 e~ %ip , 
7tt V 4 / 

1. ( HI S in 2 & - — sin 2 cos 2 e 2 ^ , 
7tt V 8 8 / 

1_ ( Hi gin 2 # - — sin 2 cos 2 e" 2 ^ , 
7tt V 8 8 J 

— I 12 sin "d cos $ sin cos 3 1? | e~ l</3 , 

7vr V 4 / 



;r~ 1 / t: ( cos 4 $ — ~r cos 2 $ + t 



J_ /3 /35 
7^ V2 VY 
1 /3 /35 



7 N 



~ \l - [ — cos 4 - — cos 2 i? + I) e~ 2iv 
7vr V 2 V 2 4 4/ 

I [S / Q 35 

— i / - I — sin i9 cos $ H sin $ cos 2 $ ) 

7vr V 2 V 4 4 

^./^ sin 3 cos tfe- 3 ^ 
7vr 4 V 2 

1/21 2 n 105 . 

— — sm 17 sin' 

7tt V 8 8 



0cos 2 <&) e~ 2i * 



— J — (^— — sin $ cos & + — sin $ cos 2 i? ) <• ' 



1 35 /3 
7^TV2 



sm 



3 1? cos tfe~ 3 ^ , 



1/3/7 21 2 n 35 4 q 

"7^V2U + t cos ^-y C0S * 



1 35 3 



7tt 8 V 2 
19 



sin 4 $ cos 0e 4j</3 . 



(66) 
(67) 
(68) 
(69) 
(70) 
(71) 

(72) 

(73) 

(74) 

(75) 
(76) 

(77) 

(78) 

(79) 

(80) 



The combinations entering Eqs. (JTSJl for ru are 

S m {<»,<p) = ReS^ + ReS™ 20 



'3/105 4 n 2.7 
cos 4 $ + 21 cos 2 ■& 



2 V 4 4 



Soi(#,¥>) = Re^ + Re^r-Re^T'^-Re^T'" 1 

105 • - .3, 



24 cos $ sin $ sin $ cos i? I cos <p 



ReS™ 22 + ReS™ 22 + ReS™ 2 '- 2 1 " - " m " 2 



^1 . 2 , 105 . 2 , 2 

— sin 17 

V 4 2 
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+ Re 

sin 2 1? cos 2 1? ) cos 2ip , 



5 , io(t?,y) 



Re5 2 3 120 + Re53 2120 



105 



12 cos $ sin $ sin i? cos 3 $ ) cos <p , 



Re5 2 3 121 +Re53 2121 



O c 212,-l 

Re 6i 



; 13 



Re S; 



212,-1 
31 



3 /35 4 63 2 7 

- (^y COS $ - — COS 1? + ~ 



/3 /35 4 n 77 2 „ 7\ 

+ \ — COS 17 COS 17 H COS 2iO . 

V 2 V 2 4 4/ 



5 12 y>) = Re 5 2 3 122 + Re S 2 3 l 22 + Re S^ 12 '" 2 + Re S 21 , 2 '- 2 

9 35 \ 

■ - sin i? cos i? H sin it 1 cos 3 $ | cos 09 

2 V 4 4 / Y 



35 3 



cos i? sin 3 19 cos 3y . 



S 20 (it» = ReS 2220 + ReS 2220 
21 „ „ 105 



sin 2 "d — sin 2 1? cos 2 19^ cos 2^ 



S 2 i(it» 



Re 5 2 3 221 + Re S 2221 - Re S™'' 1 - Re S 1 . 222 '" 1 



31 



13 



31 



/ -(—— sini9cosi9 
/ 2 V 4 

35 „ 2 n\ 35 /3 3 n 

+ — sin I? cos 17 ) cos if — — y - sin 17 cos 17 cos 3<p 



20 



S 22 (0, <p) = Re S!! 22 + Re S 3 2222 + Re 5 2 3 22 '~ 2 + Re Sff 

/3/35 4o 21 2o 7 
= "V2U COS ^"T cos ^+8 

35 /3 . 4 
8 V 2 

Further calculations are conveniently carried out in terms of the harmonic variables 
sini9cos<£>, y = sin^sinyj, and z = cosi9, which satisfy x 2 + y 2 + z 2 = 1. One has 

sin 2 $ cos 2<£> = 2x 2 + z 2 — 1 , sin 2 cos 2 $ cos 2 93 = 2a; 2 z 2 + z 4 — z 2 , 
sin cos 4<y9 = 8x 4 — 8x 2 + 8x 2 z 2 + 1 — 2z 2 + z 4 , sin $ cos 3 $ cos = xz 3 , 
cos sin 3 § cos 39? = 4^a; 3 — zx + 3z 3 a; , sin 7? cos cos cp = xz . 

Accounting for these relations, one finds 

S 00 (x, z) = (-lbz 4 + 12z 2 - l) , S 01 (x, z ) = \ ( 16xz ~ 35xz 3 ) , 

21 

S 02 ( x , z) = — (2x 2 + Qz 2 - 10x 2 z 2 - 5z 4 - 1 



1 7 l~3 

S 10 ( x , z) = - (A8xz - 105xz 3 ) , S n = -J- (z 2 + x 2 - 10x 2 z 

(48x2 - 35xz 3 - 70a; 3 ; 



S 12 (x, z) = -y r ( iS.r.: - 35./-° - TO.'"' : 
5 20 (x, z) = — ( 2.r - Or - !()./■-;- - o : 1 - 1 



( x , z) = i (48xz - 35xz 3 - 70x 3 z 



7 /3 

S 22 (x,z) = (5z 4 -8z 2 + 20x 4 - 20x 2 + 20xV + 3) 



1 



4 V 2 

We may now evaluate the combinations Rm(%, z) entering Eq. (jHHJ): 

Ro(x,z) = ^=S 02 (x,z) + ^S i(x,z) + S 00 (x,z) 

= -^L{7[2A 2 x 2 + 6(A 2 + 2)z 2 - 10A 2 xV 

- 5(A 2 + 3)z 4 - A 2 - 1] + 2A!(16xz - 35xz 3 )} , 

R x {x,z) = X 2 S 12 (x, z) + XiS n (x, z) + V6S 10 (x, z) 

= ^{t 48 ( A 2 + l ) xz ~ 35 ( A 2 + 3)xz 3 - 70A 2 x 3 z] 
+ 7Ai(z 2 + x 2 -10a; 2 z 2 )}, 

21 



R 2 (x,z) = \ 2 S 22 (x,z) + \iS 2 i(x,z) + \/6S 2Q (x,z) 

= ^{7[-20A 2 x 4 + 2(10A 2 + 3)x 2 

- 10(2A 2 + 3)x 2 z 2 + 2(4A 2 + 9)z 2 - 5(A 2 + 3)z 4 - 3A 2 
+ 2X 1 (A8xz - 35xz 3 - 70x 3 z)} . 

Substitution of Eq. (JHTJ) into Eqs. $2$-$E$ leads to the results 

Ro(x,z) = -^L{7[2A 2 x 2 + 6(A 2 + 2)z 2 - 10A 2 xV 

- 5(A 2 + 3)z 4 - A 2 - 1] + 2\ 1 (16xz - 35xz 3 )} , 

/3 

Ri{x,z) = ^={[48(A 2 + l)x^-35(A 2 + 3)x^ 3 -70A 2 x 3 2] 



+ 7\ 1 (z 2 + x 2 -10x 2 z 2 )} , 



R 2 (x,z) = ^{7[-20A 2 x 4 + 2(10A 2 + 3)s 2 - 10(2A 2 + 3)x 2 z 2 

+ 2(4A 2 + 9)z 2 - 5(A 2 + 3)z 4 - 3A 2 - 3] 
+ 2Ai(48xz - 35xz 3 - 70x 3 ^)}. 
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